
Description 

Method for determining the amplitude and phase angle of 
a measuring signal corresponding to a current or 
voltage of an electrical power supply network 

The invention relates to a method for determining the 
amplitude and phase angle of a measuring signal 
corresponding to a current or a voltage on an 
electrical power supply network by using sampled values 
of the measuring signal, a model of the measuring 
signal containing at least a sinusoidal component being 
used to calculate the model amplitude and phase 
parameters of the measuring signal with the sampled 
values by applying a recursive least -squares estimation 
method. 

A method of this type is described in an article by 
M.S. Sachdev and M. Nagpal "A recursive least squares 
error algorithm for power system relaying and 
measurement applications", IEEE Trans. * on Power 
Delivery, Vol. 6, No. 3, July 1991. In the case of this 
known method, sampled values are formed from a 
measuring signal corresponding to a current or voltage 
of an electrical power supply network and, by means of 
a linear least-squares estimation method, -by using a 
sinusoidal signal model that models the measuring 
signal, are used to form a complex vector, which 
specifies the amplitude and phase angle of the 
measuring signal. In this case, in a first step the 
- real - and - imaginary - part- of -the vec tor- -are -in each- case - 
determined on their own. Then, in a second step, by 
means of a coordinate transformation, the polar 
coordinate representation of the complex vector, that 
is to say the magnitude and phase of the vector, can be 
determined from the real and imaginary part of the 
vector . 
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In the case of the known method, it is assumed that the 
frequency of the measuring signal is known. If this is 
not so, or if the frequency changes, then a separate 
method is needed in order to determine the frequency of 
5 the measuring signal. Methods are known there, for 
example, which measure the distance of the zero 
transitions of the measuring signal and, on the basis 
of measuring this period, determine the frequency of 
the measuring signal, see, for example, E. Schrufer 
10 (Editor) : "Lexikon Me£- und Automat isierungstechnik" 
[Measurement and Automation Encyclopedia] , VDI-Verlag, 
1992, p. 204. Also known is a method for frequency 

O measurement in which the measured signal to be examined 

is in each case filtered in parallel by a high-pass 

Sj 15 filter and an all-pass filter (German patent 

H DE 4 2 11 94 6) . The frequency of the measuring signal 

can be determined via the ratio between the amplitudes 

□ of the output signals from these two filters. 

•I i 

|1J. 2 0 The invention is based on the object of specifying a 
;r[ method with which all the significant variables of the 

measuring signal can be determined simultaneously and 

quickly. 

2 5 In the case of a method of the type indicated at the 
beginning, this object is achieved, according to the 
invention, in that use is made of a model of the 
measuring signal containing the sinusoidal component in 
accordance with the relationship y = A ■ s±n(27rft + (p) , 
— - 30 ~y designating- an- i-ns-t-a-nta-neous-- value -of- the— model— of. 

the measuring signal, A the amplitude, f the frequency, 
(p the phase angle and t the time; by using this model 
of the measuring signal and by using the sampled 
values, by means of a recursive nonlinear least -squares 

35 estimation method, the model frequency parameter of the 
measuring signal is also determined by the estimation 
together with the model amplitude parameter and the 
model phase angle parameter. 
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Although it is known, from the book by H.-J. Hermann, 
"Digitale Schut ztechnik" [Digital Protection 

Engineering], 1997, pp. 110-111, to use a recursive, 
nonlinear least-squares estimation method in protection 
5 engineering, the book does not contain any reference to 
the fact that by using such an estimation method, 
sampled values of a measuring signal can be used to 
determine the amplitude, phase angle and frequency of 
the measuring signal in a single measured-value 
10 processing process. 

A significant advantage of the method according to the 
invention is, however, precisely that the sampled 
values of the measuring signal are used to determine 
15 the frequency, as well as the amplitude and phase 
angle, in a measured-value processing process, and 
therefore the amplitude, phase angle and frequency of 
the measuring signal are associated with the same point 
in time. 

20 

Furthermore, the patent DE 42 05 300 CI discloses a 
method with which the phase angle and the amplitude of 
a periodic signal can be determined via a phase-locked 
control loop (PLL (Phase-Locked Loop) ) . 

25 

The use of the model for the measuring signal 
y = A • sin(2^ft + cp) leads to good results if the 
measuring signal has a purely sinusoidal waveform. If 
there is a DC component present in the measuring 
- 3 0- signal-,- - a model of- the- -measuring, -signal - in- accordance, 
with the relationship y = A • sin(2;rft + <p) + d is 
advantageously used, the summand d modeling the DC 
component of the measuring signal . 



f 

y = A- sin 
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If the measuring signal is a signal whose frequency 
changes over time, in the case of such measuring 
signals without a DC component, a model of the 
measuring signal in accordance with the relationship 

V /=o y 
can advantageously be used, and in the case of 
measuring signals with a DC component, a model of the 
measuring signal in accordance with the relationship 

y = A-sm 2xJ](f {n t l )t + <p \ + d 

can advantageously be used, f u> designating the ith 
time derivative of the frequency and modeling a change 
in the frequency over time, and various orders of the 
time derivative of the frequency being taken into 
account by selecting the variable n. In these models, 
in an extension of the model specified above, 
y = A • sin(2;rft + cp) , the frequency 

f is replaced by the expression ^(/ ( 'V) ; if only 
a 0th order derivative is taken ,ae0 into 
account, the integral expression becomes f (0) or f. 
These expanded models make it possible to determine 
frequency changes f (l) over time, in addition to the 
variables amplitude A, phase angle q> and frequency f . 

In a further advantageous refinement of the method 
according to the invention, the values of the amplitude 
A, the phase angle (p and the frequency f determined by 
the estimation method are output as results only when 
~ the "estTma t~i~ on" error s~less t harT "a" ~sma I'l est p"e^mi"tt ed" 
estimation error. This has the advantage that, in 
particular, the values estimated at the start of the 
method and afflicted by large estimation errors are not 
output, and therefore the large 
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estimation errors cannot have any negative consequences 
for a user of the method. 

In order to explain the invention further, 
5 figure 1 shows, by using a block diagram, the progress 
of an exemplary embodiment of the method according to 
the invention, 

figure 2 shows waveforms of a measuring signal, of an 
estimated signal y s determined for the measuring signal 
10 by using a model y = A • sin(2^rft + (p) and of the 
resulting error F, 

figure 3 shows results of the method according to the 
invention in the case of a purely sinusoidal measuring 
signal when use is made of an appropriate model of the 
15 measuring signal, separated in accordance with the 
time-based waveforms of amplitude, phase angle and 
frequency, 

figure 4 shows results of the method according to the 
invention in the case of a measuring signal subject to 
2 0 an offset, that is to say containing a DC component, 
separated in accordance with the time -based waveforms 
of amplitude, phase angle and frequency, when use is 
made of a model of the measuring signal that does not 
model any DC component, and 

2 5 figure 5 shows results of the method according to the 

invention in the case of a measuring signal subject to 
an offset, that is to say containing a DC component, 
separated in accordance with the time-based waveforms 
of amplitude, phase angle, frequency and DC component, 

3 0 .when .us.e i.a _made of. a. model o.f__the measuring, signal. 

that models a DC component. 

According to figure 1, after being sampled in a sample- 
and-hold circuit 1 and after analog-digital conversion 
35 in an analog-digital converter 2, sampled values y m of 
a measuring signal u m are present on an input 3 of - a 
data processing system 4. The data processing system 4 
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contains a unit 5 for carrying out a recursive 
nonlinear least - squares estimation method. 

At the beginning of the recursive estimation method, 
5 starting values SW for the variables amplitude A, 
frequency f and phase angle cp are input into the unit 5 
and are present on the output 6 of the unit 5 as 
estimation output values . The estimation output 

values o are led to a block 7, which contains a model 

— (0) 

10 of the measuring signal. 

Q 

H In the block 7, the estimation output values are 

I"* used, takinq into account the model of the measuring 

signal, to determine a starting signal value y s0 , which 
15 is output at the output of the block 7 and led to an 
input 8 of the unit 5. 

The unit 5 uses a (first) sampled value y m of the 
measuring signal u m and the starting signal value y s0 to 
20 determine an estimation error F es tim in accordance with 
the relationship (3) recited below. This estimation 
error will lie above a permitted smallest estimation 
error. Therefore, in the unit 5 and starting from the 
estimation output values @ (Q) , in accordance with the 

25 rules for recursive nonlinear least-square estimations, 
new estimated values 0 are formed in accordance with 

the relationship (1) below, are output at the output 6 
of the unit 5 and led to block 7. 

30 In block 7, the new estimated values are used to 

determine an estimated signal value y sl , which is output 
at the output of block 7 and is led to the input 8 of 
the unit 5. The unit 5 uses the (first two) sampled 
values y m and the starting signal value y s o and the 
35 estimated signal value y si again to form an estimation 



i — = 



10 
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error F es tim in accordance with the relationship (3) 
recited below. This estimation error will generally 
also lie above the smallest permitted estimation error. 
Therefore, in the unit 5 and starting from the 
estimated values ©^ on the preceding estimation, new 

estimated values © are again determined in accordance 

~— — (2) 

with the relationship (1) . These steps are repeated 
until the estimation error F es tim lies below a smallest 
permitted estimation error. 



In order to determine whether the estimation error F es tim 
lies below a smallest permitted estimation error, use 
is made of a block 9. Via an input 10, block 9 obtains 
M from the unit 5 the sampled values y m of the measuring 

f 7j 15 signal u m , the starting signal value y s0 and the 

gj estimated signal values y si to y S k/ which are buffered 

in the unit 5. When the estimation error F est im lies 

Q 

ill below a smallest permitted estimation error, then a 

[U switching signal, which switches a switch 12 on, is 

1-2 20 output at an output 11 of the block 9. The switch 12 

111 then forwards the estimated values © present at the 

output 6 of the unit 5 to an output 13 of the data 

processing system 4. The estimated values © are then 
A — (*) 

output as the resulting values from the estimation 

2 5 method and constitute the values for amplitude, 

frequency and phase angle of the measuring signal, 
estimated with sufficient accuracy. After that, the 
estimation method continues to run as described above 
and permits continuous determination of" the "estimated" 

3 0 values © . 

— {*) 

The unit 5 determines the estimated values ©^ in 



accordance with the following relationship (1) . 

(i) 

iv„. - n\(p . ,w, t _ n // 

-(*) 
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In relationship (1) , @ represents a vector which 

V* ) 

contains the estimated values of the variables to be 
determined, that is to say, here, the amplitude A, 
frequency f and phase angle cp, after k estimation 
steps; @ (A designates a vector which results in the 

case of the estimation after k-l estimation steps. The 
matrix P(k) is a symmetrical precision matrix, as it is 
known, whose determination is presented further 
below using the relationship (2) . The function 
contains the model of the measuring signal. The 
vector Xosl contains the partial derivatives of the 
model ^^ (Jt) >®(*-n j of the measuring signal with 

respect to the parameters of the model of 

the measuring signal, that is to say the derivatives 
with respect to amplitude, phase angle and frequency. 
The variable y (k) is the kth sampled value of the 
measuring signal . 

In order to determine the precision matrix P (k) , use is 
made of a relationship (2) . 



p = — 

- (k) X 



p (k _„ — 7 y £(k-D 

" <k X + y T P< k -„Y {k) -"' 

— IK) — 



(2) 



The factor X determines the exponential weighting of 
preceding sampled values of the measuring signal . At 
the beginning of the method, a starting value is used 
for the precision matrix P (k) . 

The estimation error F es tim after the kth estimation is 
determined by means of the relationship 



(3) 
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Here, N is the number to be taken into account of the 
sampled values y m of the measuring signal u m and of the 
estimated signal values y s/ y mi is the ith sampled value 
of the measuring signal u m and y S i is the ith estimated 
5 signal value . 

The estimation error F es tim is determined by N sampled 
values y m i of the measuring signal u m and N estimated 
signal values y S i being evaluated. An appropriate 

10 magnitude for N is the quotient of the (sampling 
frequency of the sample-and-hold circuit 1) / 
(estimated frequency) , but N can also be selected to be 
larger. During the first estimation pass, the value 
used for N is the quotient of (sampling frequency of 

15 the sample-and-hold circuit 1) / (starting value SW of 
the frequency) . When, at the beginning of the 
estimation method, N sampled values y mi of the measuring 
signal u m , or N estimated signal values y s ± are not yet 
present, then only the present values are used to 

20 determine the estimation error 

estim • 

Further details relating to carrying out recursive 
nonlinear least -squares estimation methods are 
described in the documents J. Wede ; D. Werner: 

25 "EchtzeitprozeSmodelle auf der Basis von 

Parameterschatzverf ahren" [Real-time process models 
based on parameter estimation methods] , automation 
engineering series, Vol. 214, VEB Verlag Technik, 
Berlin, 1985, p.p. 30-34, pp. 44-50, 56, 57, and Hans- 

3 0 Joachim " Herrmann": " "Digi'tale Schut-ztechni-k-: -G^und-lagen., _ _ 
Software , Ausf iihrungsbeispiele" [Digital Protection 
Engineering: Principles, Software, Design Examples] , 
VDE Verlag 1997, pp. 104-113, and A. Jurisch: "Digitale 
ImpedanzmeSverf ahren auf der Basis von 

35 Identif ikationsmethoden" [Digital Impedance Measuring 
Methods Based on Identification Methods] , Dissertation, 
Technische Hochschule Zittau, 1990, Annex 4.6.3. 
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Figure 2 shows the time-based waveform of a sinusoidal 
measuring signal u m , whose amplitude, frequency and 
phase angle 
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are to be determined by using the estimation method 
according to figure 1. Also shown over time, as the 
curve y s/ are the estimated signal values y s i . . . y 3 k 
which are output at the output of block 7. Also shown 
5 is the waveform of the values of the error F which 
results from the difference between the measuring 
signal u m and the estimated signal values y si . . . y S k/ 
over time t. The waveforms shown result from an 
estimation method which is carried out by using the 
10 model y = A • sin(2/rft + <p) for the measuring signal. It 
can be seen that the waveform of the curve y s is the 
I" 4 same as the waveform of the measuring signal u m after 

i about 2 0 ms; the error F then assumes very small 

y values. The estimation method therefore permits the 

H 15 correct determination of the amplitude, the frequency 

ij and the phase angle of the measuring angle after about 

one period of the measuring signal u m . 

Q 

h* In figure 3, taking into account the waveforms shown in 

LH 20 figure 2, the estimated values © ... © linked with 

i u ~ — (i) — (*) 

[jj the estimated signal values y s i . . . y S k are plotted 

against the time t, separated in accordance with the 
variables A, co and 27if and (p that determine them. In an 
upper graph, figure 3 shows the waveform of the 

25 amplitude A determined by the estimation method 
according to figure 1, in a central graph the waveform 
of the variable co = 27tf that is proportional to the 
frequency f, and in a lower graph the waveform of the 
phase angle cp, in each case plotted against time. 

3 0 According to figure" 2 , "the result that" the time-based 
waveforms have reached a steady state at their final 
value after about 20 ms with the smallest permitted 
estimation error, and indicate the result of the 
estimation method . 

35 

Figure 4 shows, in a manner of representation which is 
the same as figure 3, the time-based waveforms 
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of amplitude A, frequency f and phase angle cp 
determined in a further estimation by using the model 
of the measuring signal y = A • sin(2^rt + <p) if use is 
made of sampled values of a measuring signal which 
contains a DC component . The waveforms of frequency and 
phase angle also contain large time-based fluctuations 
for times greater than 20 ms . These fluctuations point 
to an erroneous estimation and are a reference to the 
fact that the selected model is not beneficial for a 
measuring signal containing a DC component. 

Figure 5 shows the time-based waveforms of amplitude A, 
frequency f, phase angle cp and DC component d 
determined in an estimation using the model 
y = A • sin(2/rft + cp) + d for the measuring signal, if 
use is made of sampled values of a measuring signal 
which contains a DC component. By using the model just 
mentioned of the measuring signal, fluctuations in the 
time-based waveforms of frequency and phase angle are 
avoided; consequently, an estimation which is better 
with regard to the estimation error is carried out. 



